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Abstract 

We propose a model and an estimation technique to distinguish sys¬ 
temic risk and contagion in credit risk. The main idea is to assume, for 
a set of d obligors, a set of d idiosyncratic shocks and a shock that trig¬ 
gers the default of all them. All shocks are assumed to be linked by a 
dependence relationship, that in this paper is assumed to be exchange¬ 
able and Archimedean. This approach is able to encompass both systemic 
risk and contagion, with the Marshall-Olkin pure systemic risk model and 
the Archimedean contagion model as extreme cases. Moreover, we show 
that assuming an affine structure for the intensities of idiosyncratic and 
systemic shocks and a Gumbel copula, the approach delivers a complete 
multivariate distribution with exponential marginal distributions. The 
model can be estimated by applying a moment matching procedure to the 
bivariate marginals. We also provide an easy visual check of the good spec¬ 
ification of the model. The model is applied to a selected sample of banks 
for 8 European countries, assuming a common shock for every country. 
The model is found to be well specified for 4 of the 8 countries. We also 
provide the theoretical extension of the model to the non-exchangeable 
case and we suggest possible avenues of research for the estimation. 

Keywords: Credit risk, Systemic risk, Contagion, Copula functions, 
Marshall-Olkin distribution, Financial crisis 


1 Introduction 

The purpose of this paper is to draw a line between systemic risk 
and contagion, and to design a method to measure the relative con¬ 
tribution of contagion and systemic risk to the dependence structure 
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of a set of credit exposures. The statistical problem of disentangling 
systemic risk and contagion is of utmost relevance for economic pol¬ 
icy. In several problems, such as pollution regulation or banking, 
recognizing systemic risk, as an event independent of the agents, or 
contagion, that is a system wide event triggered by one of them, is 
a major discriminant factor to decide whether the effects should be 
charged to the community at large or to the individual agents. 

The problem is involved for two reasons. The first, already ad¬ 
dressed in many studies on the subject, is that the systemic risk 
factor is not observed and for credit risk applications we are only 
able to extract the marginal survival probabilities from the market. 
The second reason, that is the subject of this paper, is whether such 
dependence is explained by the presence of a systemic risk factor 
only, or of some infectious elements in the system. 

To make the problem clear, assume that we are allowed to observe 
the systemic risk factor, and we are able to appraise the probability 
of a systemic crisis. In this situation, the question would natu¬ 
rally arise whether the systemic shock is independent of the other 
events triggering the default of each component of the set. Answer¬ 
ing this question on practical grounds would be obviously easy in 
this setting, and the dependence between idiosyncratic and systemic 
triggers of default could be estimated in the usual way, e.g. using 
copulas. 

Notice, however, that even if dependence with the systemic risk fac¬ 
tor were observed, a problem of interpretation of these results on 
theoretical grounds would arise concerning how it would affect the 
observed dependence structure of the components in the system. 
In fact, the presence of a systemic factor is sufficient to induce de¬ 
pendence among the components of the system and between each 
component and the systemic event. In other words, this depen¬ 
dence shows up even if the systemic shock is independent of the id¬ 
iosyncratic ones. Intuitively, if the idiosyncratic default drivers were 
linked to the systemic risk trigger by a dependence relationship, the 
degree of dependence in the system would be even stronger. 

In this paper we propose a model to represent these two sources 
of dependence in a tractable way. The idea of our model is very 
simple. Given a cluster of d obligors, we assume that the system 
is subject to a set of d + 1 shocks, one of which is common to all 
the components and leads to simultaneous default of all the oblig- 
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ors in the cluster, while the others are responsible for the default 
of each component. Assume further that the times of these shocks 
are linked by a copula function of dimension d + 1. It is immediate 
to see that this model encompasses the two extreme cases of pure 
systemic risk and pure contagion. Namely, in the case of a product 
copula for the d + 1 occurrence times of the shocks, we obtain a 
model with bivariate Marshall-Olkin marginal distributions, repre¬ 
senting the pure systemic risk model. The opposite case arises when 
the systemic shock has zero probability, so that we have a standard 
survival copula model corresponding to pure contagion. Allowing 
for positive probability of a systemic crisis and for dependence be¬ 
tween this risk and the idiosyncratic credit drivers would then allow 
to design models that represent both systemic risk and contagion. 
In these models, the task of disentangling the two is a relevant ques¬ 
tion. 

In its simplest version, our paper assumes the standard restriction 
of the choice of exchangeable copulas for the credit risk drivers. 
This means that each idiosyncratic factor is assumed to be linked 
by the same dependence structure to the systemic risk factor and to 
the other idiosyncratic ones. We also show that further restrictions 
change the copula model in a multivariate distribution model with 
exponential margins. If the model may seem restrictive, on practi¬ 
cal grounds we provide a methodology to verify if the assumption 
is borne out by the data. Moreover, on theoretical grounds, we will 
also provide the theoretical development for the non-exchangeable 
version of model. 

The plan of the paper is as follows. After reviewing the relevant lit¬ 
erature, in Section [2] we motivate and describe in full generality our 
credit risk model for a basket of issuers, we discuss its main prop¬ 
erties and the restrictions that change the copula model in a mul¬ 
tivariate distribution with exponential marginals. In Section [3] we 
discuss the theoretical features of the extension to non-exchangeable 
dependence of the credit risk drivers. Finally, in Section [4] we illus¬ 
trate our application to the banking system of a set of countries of 
the Euro area. In Section [5] we report conclusions and a discussion 
of the main issues left for future research. 
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1.1 Related literature 


Our paper is related to a large literature on the measurement of sys¬ 
temic risk and contagion, even though to the best of our knowledge 
it is the first attempt to disentangle the two. Leaving aside any hope 
of being exhaustive, we may provide a taxonomy of the main con¬ 
tributions according to the structure of models and the data used. 
As for the methodology involved, a first class of models are based 
on the application of Granger causality, and related concepts, to the 
prices of financial assets (Billio et ah,2012). A second set of models 
is based on the network representation of the relationships among 
financial institutions (Diebold and Yilmaz, 2011). A third approach 
is based on the theory of risk measures applied to systemic risk and 
contagion. Models in this class are based on the measurement of 
expected losses conditional on an extreme scenario of some systemic 
risk factor. The technique is the same as expected shortfall, with the 
difference of conditioning with respect to a systemic variable. These 
measures are called Marginal Expected Shortfall, MES (Acharya et 
al., 2010), and Co VaR (Adrian and Brunnermeier, 2011). Cherubini 
and Mulinacci (2014) give conditions to ensure that coherence re¬ 
quirements be met, and propose examples of measures in this class 
based on copula functions. 

Coming to the kind of data that are used in the empirical analysis, 
we may distinguish between applications that rely on the analysis of 
market prices, and those that use flows and balance sheet data. The 
first choice use equity stock prices (Billio at ah, 2012) or volatili¬ 
ties (Diebold and Yilmaz, 2011), credit spreads of bonds and credit 
derivatives (Baglioni and Cherubini, 2013). With this choice, the fo¬ 
cus is on measurement of the effects of systemic risk and contagion, 
in terms of future cash flows and the default probability that are 
implied in market quotes. The second choice exploits flows among 
the financial intermediaries and the focus is more on the means that 
explain propagation of the shocks through the financial intermedia¬ 
tion system. Here the analysis is focussed on flows in the interbank 
market (Bonaldi, Hortacsy and Kastl, 2013) or on several layers 
representing other markets (Bargigli et ah, 2013), or else on balance 
sheet indexes such as leverage (Brownlees and Engle, 2010). All 
these proxies are used as measures of the strength of contagion in 
the system. 

Our paper uses the default probability extracted from CDS and their 
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dependence structure in order to recognize how much of this depen¬ 
dence is due to relationships among the components of the system, 
as in network based models, and how much of the co-movement is 
due to the presence of a systemic risk factor, as in systemic risk 
models. Moreover, to the best of our knowledge, this is the first at¬ 
tempt to include a dependence structure between each component 
and the systemic shock, although being the same dependence struc¬ 
ture in both cases. 


2 The model 

Here we introduce the motivation of the model, and its basic set¬ 
ting. The idea is that in a system of d components, the lifetime 
of each of them can come to an end either for idiosyncratic or sys¬ 
temic shocks, as in a standard Marshall-Olkin setting. Differently 
from that model, in which all shocks are assumed to be independent, 
here the idiosyncratic components are infectious. Idiosyncratic de¬ 
faults can be associated, and they may also represent triggers of the 
systemic shock, leading to default of the whole system. 

Technically, let (0, T ', P) be a probability space with a d + 1-vector, 
(X 0 , Xi ,..., X d ) whose components have [0,+oo) as support. X 0 
denotes the arrival time of the systemic shock and (Xi ,..., X d ) are 
those of the idiosyncratic ones. We assume that the joint survival 
dependence structure is represented by a strict Archimedean copula, 
that is 

F(x o, xi ,..., x d ) = (^~ 1 {F 0 (x 0 )) H-f ^~\F d )) 

for (x 0 ,..., Xd) G [0, +oo) d+1 , where F, (that is assumed to be con¬ 
tinuous and strictly decreasing) is the marginal survival function of 
Xi and ^ is the generator of a strict d+ 1-dimensional Archimedean 
copula. We recall that ^ is the generator of a d + 1-Archimedean 
copula if and only if : [0,+oo) —>■ [0,1] is d + 1-monotone on 
[0, +oo) that is 

• it is differentiable on (0, +oo) up to order d— 1 and the deriva¬ 
tives satisfy (—1 ) k ^ k \x) > 0 for k = 0,1,..., d— 1 and 
x G (0, + 00 ) 

• (—1 ) d ~ 1 'i/j^ d ~ 1 '> is non-increasing and convex in (0, + 00 ). 
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(see McNeil and Neslehova, 2009, for more details on multidimen¬ 
sional Archimedean copulas). 

Since we restrict ourselves to the strict case, we assume 0(x) > 0 
for all x G [0, +oo). Let us define 

r fc = min{A' 0 , X k }, k = 1,..., d. 


This is the standard Marshall-Olkin setting in which the only com¬ 
mon shock taken into account is the one affecting all the components 
in the set. Of course, other specifications are possible, including 
models with more than one common shock, affecting selected sub¬ 
sets of the components (see, for all, Durante, Hofert and Scherer, 
2010). The observed default times r k represent the first arrival time 
between a common (systemic) shock affecting all the system and the 
idiosyncratic shocks. We then add an Archimedean type of depen¬ 
dence among the arrival times of the shocks, in order to represent 
contagion. 

The joint survival function of the random vector r — (ti, ..., rf) 
can be easily recovered 

Fr(h, ... ,£ d ) = 0 \ ^ _1 (F 0 (max{4})) T V] 0 _1 (P fc (4)) ) (1) 

\ l<fc<d z —* / 

\ k =1 / 


for ti,..., td G [0, +oo) d , while the marginal survival functions are 


F Tk (t) = 0 (0 1 (F 0 (t )) + 0 \F k (t))) = 0 (H 0ik (t )), t G [0, Too) 

( 2 ) 


where H 0tk (x) =0 1 (F 0 (x)) + 0 1 (F k (x)). 

It is also easy to extract the copula function of the observed default 
times 


Proposition 2.1. The survival copula C of the vector of default 
times t is, for u G [0, l] d , 


i=i 


0 


-1 


>0 + 

k=l ,k^j 


-1 


K) !t( u ) (3) 


where D k (x) = 0 1 o F k o H 0 l(x ) and 


Aj = 



e [0, l] d : ™x{ H o,i ° 0 ^Wi)} = H o,j 


° 0 1 ( u j) 
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with the convention that if u satisfies the required condition for more 
than one index j, it is assumed to belong to the Aj with the smallest 
index j. 

Proof. See Appendix [HI □ 


2.1 A multivariate distribution with contagion and expo¬ 
nential marginals 

In practical applications it is common to represent and calibrate 
default times by exponential distributions 

F Tk {x) = exp(—/i fc x), (4) 

where /ifc denotes the intensity parameter. We now discuss which 
restrictions can be imposed on the model in order to transform the 
copula model above in a multivariate distribution with exponential 
marginals. 

Starting from the copula model illustrated, constructing such 
multivariate distribution would imply a data generating process such 
that: i) the distortion functions Dk(x) are linear; ii) the dependence 
is represented by a Gurnbel copula. 


2.1.1 Linear distortion 


A possible assumption about functions in the spirit of 

the paper by Muliere and Scarsini (1987), is that they are all pro¬ 
portional to the same function K(x): that is, ^” 1 (F i (x)) = A iK(x) 
for A* > 0, for i = 0,1,..., d. This is equivalent to Dfix) = (1 — ctj)x 
where 


An 


OLi = 


6 [ 0 , 1 ) 


A* + Ao 

and the obtained copula is independent of K. 

In the more specific case in which if is completely monotone (that 
is is the Laplace transform of some positive random variable), 
we recover the Scale-Mixture of Marshall-Olkin distributions and 
copula models (SMMO) studied in Li (2009). The exchangeable 
case of SMMO model is studied in Mai and Scherer (2013) where it 
is applied to the pricing of CDOs. 
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2.1.2 The Gumbel case 


A further restriction to yield marginal exponential distributions is 
to consider the case in which ^ is the Gumbel generator, that is 


'ijj(x) = e xV , 6 > 1. Now, equations (□]), (EJ) and (j3]) take the form 


F r {t i,..., t d ) = exp ^ - ( A 0 K l max{i,} j + ^ A k K(t k 


k =1 


F Tk {t) = exp ( -(A 0 + Afc) e K 6 (t) 


(5) 


e( u ) = E 


exp 


(-In Uj) e + ^ (1 - a k )(-\nu k y 


Wu) 


3 = 1 


k=l,k^j 


Notice that, in this case, i/; is the Laplace transform of an ^-stable 
distributed random variable and so it represents a specification of 
the SMMO model of Li (2009). 


Notice that setting K{t) = t e in (|5]) yields exponential marginals 
as required 

A l k — (Ao + A/.) 9 (6) 

where //k is the intensity in equation (J4}. 


2.2 Properties of the model 

The main feature of our model, right from the most general setting, 
is to increase the degree of dependence among the default times, 
both with respect to the standard Archimedean copula without any 
systemic risk factor and the Marshall-Olkin copula in which the 
systemic risk factor is independent of the others. 

The dependence structure of the model encompasses both the 
sensitivity of the default times to the systemic shock, and the de¬ 
pendence among the shocks, represented by Archimedean copulas. 
Both these elements interact to determine the dependence among 
default times. 
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In the general setting, the Kendall’s tau measuring the de¬ 
pendence of the pair of default times ( Ti,r k ) can be written as 

pOO 

Tj, k = + 4 / ■ T(x)dx 

Jo 

where denotes the Archimedean Kendall’s tau corresponding to 
the generator i/j and 

T(x) = ° F 0 o (t/A 1 oF 0 + V* -1 ° Fj + ir 1 o F k ) 1 (x) 

where we refer the reader to Mulinacci (2014) for the derivation. No¬ 
tice that if we are interested in representing the dependence struc¬ 
ture between the systemic shock and default times, we have that 
the Kendall’s tau r^o of the pair ( Tj,X 0 ) is 

/ OO / v 

(^'(x)) 2 -(V’ -1 O F 0 O (if; -1 O F 0 + i/j~ 1 O Fj) 1 (x)j dx 

The first term is simply the Kendall’s tau of the Archimedean cop¬ 
ula used in the analysis, while the other term, that is more complex, 
involves both the generator of the Archimedean copula and the rel¬ 
ative relevance of systemic and idiosyncratic shocks. 

In the multivariate distribution arising with linear distortions and 
the Gurnbel copula in the model, these relationships simplify sub¬ 
stantially. In fact, let Cj tk (u,v) be the general marginal 2-copula, 

Cj,k{u,v) = exp |- [(-lnu) 0 + (1 - a k )(- lnx) 0 ] 8 j lp^-i (u)>a k ^(v)}+ 
+ exp | — [(1 - aj)(-lnu) 9 + (-Inn) 0 ] 9 j l{ aj ^(u)<a k ^(v)} 


Since this family of copulas represents a particular specification of 
the Archimax copulas of Caperaa et al. (2000) and of the Archimedean- 
based Marshall-Olkin copulas of Mulinacci (2014), its Kendall’s tau 
is known to be 


6-1 

T j,k ~ — 7 — 


t mo 

T j,k 

9 


( 7 ) 


where 


MO 

T j,k 


OLjCX-k 

Oij OLfc — QLjOLiz 


is the Kendall’s tau of the Marshall-Olkin copula. 
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Now, the dependence between each default time and the time of 
a systemic shock is linear 

6 — 1 a- 

T °'i = — + T 

This relationship will be used in our estimation strategy in order to 
verify the specification of the model. 

2.3 Estimation strategy 

In the estimation of the model we assume to observe a panel set of 
data representing marginal default intensities of k — 1,2,... ,d 

components, for {ti,t 2 , ■ ■ ■ ,t m } dates. Our task is to estimate the 
set of parameters, representing the sensitivity of each obligor to 
the systemic shock, and the parameter 6, that measures the degree 
of contagion in the system. We also would like to make a check of 
the specification of the model. 

Since the main feature of our approach is to identify the weight of 
the sensitivity to the systemic shock and of the degree of contagion 
in the dependence structure of default times, a natural estimation 
strategy would be a moment based approach, which resembles the 
calibration procedure proposed by Genest and Rivest (1993). In 
particular, our model specification based on linear distortions and 
Gumbel dependence makes a procedure based on Kendall’s tau cal¬ 
ibration very easy. 

Since the model is built to be fully characterized by the bivariate 
marginals, the estimation is naturally performed by calibrating the 
bivariate Kendall’s tau statistics of the system. For each cluster that 
we expect to be part of the same exchangeable system, consisting of 
d units, we calibrate the set of d + 1 parameters of our model. 

Formally, we first estimate the Kendall’s tau statistics of all the 
pairs of the sample, and then estimate the set of parameters © = 
{«!, a 2 , ■. ■, a dl 9} by solving 


d— 1 d 

© = argmin .IE dtst(fij, Otj, $)) 

{ai,a 2 ,...,a d ,0} i=1 j=i+1 

where dist(x,y ) is a suitable distance measure, Tjj(cp, «j, 0) is the 
theoretical Kendall’s tail based on estimates, and fis the corre¬ 
sponding empirical Kendall’s tau statistics. As for the parameters 
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set, cq represents the sensitivity of component * to a systemic shock, 
and 9 represents the contagion parameter, that is assumed to be the 
same across all pairs. 

The structure of the model also provides an easy procedure to 
check whether the specification of the model provides a good fit to 
the data. The idea is that if the model is well specified we could 
use it to estimate the intensity of the systemic shock from market 
data, and use that information to directly verify the specification of 
the model. The property of exponentially distributed marginals of 
the Gurnbel specification is particularly useful in this case. Given 
a panel of m observations of d intensities, we can estimate a time 
series of m intensities of the systemic shock. Using equation (JHJ) and 
the definition of cq, it is straightforward to compute 


Ao (U) 


Efc=ihfc(^) 

Z-jk=l ah 


(9) 


where Ao(tj) denotes the estimate of the systemic shock intensity at 
time t % . 

A straightforward visual check of the specification of the model 
would then be to estimate the Kendall’s tau value between the ar¬ 
rival time of a systemic shock and marginal default times. If the 
model is well specified, the Kendall’s tau values should be aligned 
on the straight line described by equation (JSJ) . 

In that case, the procedure also provides a new series representing 
the implied intensity of the systemic shock, that may be usefully 
applied for further investigation of the cluster and of the system as 
a whole. As an example, one could verify whether other elements of 
the system, originally not associated to that cluster, actually have 
the same dependence with the systemic shock as the other elements 
of the cluster. As a second example, one could use the estimated 
systemic shock intensities of different clusters to check the degree of 
association across clusters. 

Our application, that is meant to illustrate this estimation proce¬ 
dure, will be focussed on a set of European banks. We will assume 
that the banks of the same country constitute a cluster, and we will 
verify in which case this assumption is borne out by the data. 
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3 An extension to hierarchical Archimedean risk 
factors 

In this section we will consider a possible extension of the model 
with exchangeable dependence structure presented above. Clearly, 
any d + 1-dimensional copula can be considered in place of the 
Archimedean one and the same construction implemented. Among 
the possible reasonable choices, vine- Archimedean copulas and hier¬ 
archical copulas (HAC) could be considered as natural non-exchangeable 
extensions. 

In this paper we will consider d+1- dimensional HAC copulas. These 
are obtained through the composition of simple Archimedean copu¬ 
las: such composition is recursively applied using different segmen¬ 
tations of the random variables involved. Starting from the initial 
variables u±,, u d +i, these are grouped in l\ copulas ..., Ci^. 
Then, these copulas are grouped in l 2 copulas C 2 , 1 , • • •, C^, and up 
to the last level where we have just one copula. In order to ensure 
that the so obtained HAC copula is indeed a copula, the generators 
ipi.j of the copulas involved have to be completely monotone and 
the same must hold for their compositions Vd+i ,j ° ipi.k whenever C t) k 
is an argument of C l+ \, r When the generators ipij are in the same 
parametrized family, the described procedure yields a copula if inner 
copulas have a parameter higher than the outer ones: in this paper 
we will consider generators belonging to the same family (see Savu 
and Trede 2008 and McNeil 2008 amog the others as references on 
this topic). 

In the fully nested case we have 

C(u) = C d (...C 3 (C 2 (Ci (it i,u 2 ),u 3 ) ,u 4 ),,u d+ 1 ) . 

If the probability distribution of the systemic shock Xo corresponds 
to u i, then, the idiosyncratic risks A, : , i > 1, can be decreasingly 
ordered with respect to the dependence to X 0 being 

C x 0 , Xi (u,v) = C'j_i(u, v). 

If, instead, the probability of X 0 corresponds to u d+ 1 , then 

C Xo ,Xi(u,v) = C d (u,v). 

and the dependence structure between each idiosyncratic risk and 
the systemic one is the same for all the idiosyncratic triggers. 
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In the intermediate case in which the probability of X 0 corresponds 
to Uj for some j — 2,..., d, we have that 

C Xi ,X 0 (u,v) = v) 

for those Xi whose probabilities correspond to those Ui with i < j, 
and 

C Xo ,Xi(u,v) = Ci_i(u,v) 

for the probabilities of those X, L correspond to Ui with i > j. 

Of course, under other hierarchical configurations, completely dif¬ 
ferent relationships among the systemic and the idiosyncratic risks 
can be modelled. For example if 

(T ( U) C (C htl {U\ , . . . , Uj — 1) , Ch,2 (Xj 5 U j+ 1 ) 1 ' ' > 'U'd+l )) 

where C) lA and C) h2 are again HAC copulas, and X 0 corresponds to 
Uj, we have that 

C Xi ,x 0 (u,v) = C(u,v ) 

for all probabilities Xi that correspond to those Ui with i < j and 

C x 0 ,Xi(u,v) = C ht2 (u,v) 

for Xi probabilities that correspond to w, with i > j. Hence, in the 
first case, the dependence structure between Xi and X 0 is constant 
and weaker than that in the second case where however it varies 
according to the structure of Ch,2- 

Notice that, however, whatever is the case, the dependence struc¬ 
ture between A'o and Xi is always Archimedean, exactly as in the 
exchangeable case investigated in Section [2j As a consequence, the 
formulas there presented for the Kendall’s tau between the sys¬ 
temic shock and every default time continue to hold. In particular, 
if: i) all the copulas involved in the hierarchical construction are 
of Gumbel type and ii) for every idiosyncratic shock arrival time 
X t there exists a function AT* such that F 0 (t) = ^(XojK^t)) and 
Fi(t) = ^(A iKi(t)), then (jHJ) applies. Moreover, for all those default 
times Tj such that the corresponding idiosyncratic shock arrival time 
Xi has a dependence relationship with the systemic shock one X 0 
expressed by the same Gumbel copula , the pairs (ay,T 0 y) must lie 
on the same straight line (151). 
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Notice, then, between the fully exchangeable system, and the 
fully non-exchangeable one, we can identify an intermediate case 
in which the exchangeability concept is only applied to the bivari¬ 
ate relationships between the systemic shock arrival times and the 
idiosyncratic shocks, whatever the dependence among the idiosyn¬ 
cratic shocks could be. 


3.1 Dependence structure of observed default times 

However, since the statistical procedure presented in Section 12.31 is 
based on the estimation of the pairwise dependence structure of the 
default times Tj, we will now compute the Kendall’s function and 
the Kendall’s tau of any pair of default times. 

Clearly, the shocks involved are the systemic one and the two 
idiosyncratic ones that correspond to the default times we are con¬ 
sidering. Formally, let X i} Xj, X k be the three shocks arrival times 
we are considering. Whatever the hierarchical structure is, their 
joint survival distribution is of type 

F(xi,Xj,x k ) = C ^ (Cp e (FfoJjFjfe)) ,F k (x k )) 

where CW and C^ g are bivariate Archimedean copula functions with 
generators ^ and i pg. 

Here below we drop the notation according to which the systemic 
shock arrival time is denoted X 0l so that we can move it in different 
places of the hyrarchical structure. In particular, it is sufficient to 
study the two cases: the systemic shock is represented by X\ and the 
case in which it is represented by X k . For the sake of simplicity, we 
will assume that all marginal survival distributions are differentiable 
when needed. 

3.1.1 Xj is the arrival time of the systemic shock 

Assume X t be the systemic shock’s arrival time and 
Tj = min (Xi,Xj), r k = min(Xj,X fc ) 
be the considered default times. Then 

F Tj , Tk (tj, t k ) = Vty fe 1 ° {ipg 1 ° Fi(max(tj,t k )) + ipg 1 o Fjitj)) + ipj 1 o F k (t k )) , 
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F Tj (t) = ipg (i’ e 1 o Fi(t) +i> e 1 o Fj(t )) = i> 9 ° H 0jj (t) 

and 

Fr k {t) = ^ (V^ 1 O Fi(t) + 00 1 o F k (t)) =^o H 0tk (t) 

where H 0jj (t) = ipg 1 o F^t) + i/jq 1 o Fj(£) and H 0tk (t) = V^ 1 o F0£) + 
o F k (t). Hence, thanks to Sklar’s Theorem, from 

tj = H o] ° ^e 1 ( y u i) and tk = H o,k ° 00 1 (“fc) 
we get that the associated survival copula is 

^Tj,T k ( U jj U k) ~ 

= 00 ( 00 1 ° 0 e (00 1 ° Fi(max(Hoj o H~l o ip^(u k )))+ 

+00 1 o Fj o H~j o ipg'iuj)) + 00 1 oF k o H~l o ^\u k )) . 

Set 


Dij = ip e 1 oF i oH 0 j, D ik = ip 0 1 oF i oH^ k , D j{ = ip e 1 oF j oH 0 ^ D ki = ip/oF^H^. 
Then 

Aj ,T k ( U j ) u k ) = 

= 00 ( 00 1 0 0 e (max(Ai ° 0 e " 1 (%), D ik o ^(it*,)) + % o ^ 1 (%)) + T> fci o 

_ f 00 (00 X (%) + D k i ° ^(ufc)) , u k > h(uj) 

I 00 (00 1 0 0 e (A fc 0 00 1 («fc) + Dji o 0 e " 1 (n J )) + D ki o ^ («*,)) , u k < h(uj) 

( 10 ) 

where 

/i(x) = Vty o A 0 10 A? 0 0 e" 1 (^)- 


Restriction on the distribution of X i: Xj,X k 
Assume that there exist two functions Ji and K such that 
07 1 o F*(t) = Ai^(t), 0,' 1 o £,-(*) = 

and 

00 1 o F 0 £) = Aj/T(t),' 00 1 o F fc (t) = A fe AT(£) 

which implies that 

AT(0 = 0 00 (AiAT(0) • (11) 
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Now, setting /i, ;j = A* + A j and p, ik = Xi + X k , 

= fiijK{t ) and H 0jk (t) = n ik K(t) 

and 

Dij(x) = —x,D ik {x) = —x,Dji(x) = —x,D ki (x) = —x, 
ftij ftik ftij ftik 


from which 

F Tj ( t ) = ip 9 

and 

^Tj,T k (Uj,U k ) = 

'.-1 


and F Tk (t) = ip# (/x ik K(t )) 


\ j — i / \ A j 


A 


= % ( V’ 0 1 ° I max I —ip e 1 (u j ), —V’/K) I + —ip e 1 (u j ) I + —V’/K) 




ftik 


Hij 


ftik 


( 12 ) 


Remark 3.1. The Gumbel case 

1 A 

Assume ipg(x) = e _a;F and ip$(x) = e ~ x ? unth 6 > <p > 1 . T/ien 
ipf 1 o ipg(x) =i® and / f7H) writes 

^Tj,T k { U j} U k) ~ 


= exp 


max [ — (-log(uj)) 0 ,— (-log(u fc ))^ J + — (-log(itj )) 0 
ftij ftik I ftij 


+ —(-log(u fc ))* 

ftik 


+ 


Necessarily, by / fill) . X,K = A f K* and an admissible choice is 


e 


A i = A f and K = K% 


In particular, if K(t) = t^ and K(t ) = t e we recover exponential 
marginal distributions, that is 


F t .( t) = e~W and F Tk (t) = e“W. 
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The Kendall’s function and Kendall’s tau 


Theorem 3.1. If p = 0^ 1 o fj e , let (see ( flflj) ) 

C(u, v ) = 

A + Dki o ^f l (v)) , u > h(u) 

4ict> (P (Afc ° 00 1 A) + Dji o 00 1 A)) + D ki o 00 1 A)) , v < h(u ) 

where 

h{x) = A ° D" 1 ° Aj ° AA A)- 

We /lave that the corresponding Kendall’s function JC(t) = P(C(u, i>) < 
t) and Kendall’s tau are, respectively, 

rDi k oG- 1 <y<i>- 1 (t) 

D ki o^-\t) 


m = A 




and 


P'°P A ~ D ki o D ik (z))dz 
(13) 


r = 1 + 4 / ^ o 0 0 (t) ■ D ki o 0 0 (t)dt 
Jo 

rDikoG- 1 ^- 1 ^) 


Afc°V , 0 1 (*) 


P ° P (% A - D ki ° Afc (■ z))dzdt 


-4 / A 

Jo 

where 

G(x) = iff 1 o A ° Aj 1 ° AfcA) + A; A)- 

Proof. See Appendix [61 


□ 


Remark 3.2. The Gumbel case 

In the setting of Remark \3.1[ we get 


t 


£(t) = t--(-logt) 1 *■ 

0 

and 




l-^ik 0 I ^ 


-i-(-log t)* 
P%k 


logt)+-^z\ dz 
X / 


T = l+ ~V~e /w 

pik 0 V Jo 



A G -i( ( _ logt} 0) 


(-logt)^ 


log A 


1—S 


At; ) dz | dt 

X 


where G(x) — ( — ) s« + —x. 

' fJ'ik ) (J'ik 
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■6-1 








3.1.2 X k is the arrival time of the systemic shock 

Here we assume that X k is the arrival time of the shock and 
Ti = min(X i,Xk) and Tj = min(Xj, X k ) 
the considered default times. Then 

Fn^iUXj) = 00 [% l ° 00 (00 1 (^(0)) + 00 1 (A(0))) + 00 1 (A(max(0,0))] 

and 

F Ti (ti) = 00 ° ) and F Tj .( 0) = 00 o H 0J (tj). 

where 


tf<M = 00 1 ° A + 00 1 ° A and H 0>j = 00 1 o Fj + 00 1 o F k . 

!f P = 00 1 ° 00, 

A,r,(0,0) = 00 [p (p _1 o 00 1 o ^(0) + p' 1 o 00 1 o 0,(0)) + 00 1 ° Ffc(max(0,0))] 

and, applying Sklar’s Theorem, we recover the associated survival 
copula is 


CruTjiUi, Uj) = 00 [p (p 1 O 0 0 1 O Fi o O 0 0 1 (m0 + P 0^000 # 0 j 

+00 1 ° 0fe(max(0f0 l 1 O 00 1 (u i ), Adj o 00 1 (%)))] . 

Set D ik = 00 ' o0 ? oP/,0 1 , D jk = 00 :1 oFjoH-j, D ki = ip~ l oF k oF[~} 
and D k j = 00 1 o F fc o It follows 


°00 1 K)) + 


C Ti ,Tj(uii Uj) — 

00 [p (p _1 ° Afc ° 00 1 («i) + P” 1 ° Djk ° 00 1 ('Uj)) + max(D fe j o 00 1 (w0, D kj o 00 1 (u J -))] = 
f 0* [p (p _1 ° Aa, ° 00 1 («i) + p _1 ° Afc ° 00 1 (%)) + Aw ° 00 1 ( u »)], u j > 

\ 00 [p (p _1 O Djfe O 00 1 (u0 + p _1 O D ifc O 00 1 (u i )) + D k j O 00 1 (m j -)] , Uj < h(ui ) 

(14) 

where 

h(x) = 00 o D-j o o 00 1 (x). 


Restriction on the distribution of Xj, Xj, X k 

Assume there exists a function K such that 00 1 o F v (x) = X v K(x) 
for v = i,j, k, and set jjL ik = Xi + X k and Hj k = Xj + X k . It follows that 
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D ik (x) = -j^x, D jk (x ) = -j^x, D ki (x ) = Fl. x and D kj {x) 
and the marginal survival distributions can be written as 

FtM = Vty {UskKit )), s = i,j 
while the associated survival copula as 

F'Ti,Tj( U il U j) = 


= —x 

f-Ljk 


A p(p 1 (— tf’/Ai)) +P 1 (—■ V’/Ai))) +max(— • 0/(ui),— ‘tf’/Aj)) 
V \l^ik / \f^jk / / ftik f^jk 


Remark 3.3. The Gumbel case 

1 i 

Assume ipo(x) = e ~ xF and ip<j,(x) = e~ x<t> , with 9 > <f> > 1. Then 
p(x) = x% and 


CVj ,Tj A/ j bi) 



V*’ 

l^jk 


(-iogK)) 0 ] + 



(-log ko: 

l^jk 

while, if K(t) = tP, we get exponential marginal distributions 

1 1 
F T .(t) = e-*** and F T .(t ) = e"Ad 


The Kendall’s function and Kendall’s tau 


Theorem 3.2. Let (see\Tf\) 


C(u, v ) = 

f i>cj> [p (p _1 ° Afc ° A) + P _1 o D ifc O i’f l {v)) + D ki O ^ A)] , V > h{u ) 

~ \ A> A (p _1 ° A* ° V ^ 1 A) + p _1 ° Afc ° Vv 1 A)) + D kj ° A)] , v < h(u ) 

where 

h(x) = ^° Djy o o 'i/jf 1 (x). 
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We have that the Kendall’s function lC(t) = P(C(-u, v) < t) and the 
Kerndall’s tau respectively are 


K(t) = t + if'^oif^ft) 


rp 1 ojp r/> 1 oF j oG 1 o ^(t) 


I p- 1 oD jk oip (j> 1 (t) 


p' ° P 1 {V’/W - ° F k ° Fj 1 O ^ O p(z)} cfe+ 


Ap 1 {% 1 (*) -% 1 °F k oF i 1 o^o p{z)} dz+ 


J p~ 1 oD ik oi/j it) 1 (i) 

~(D tj + D ki ) a ft‘(t) + 2ft 1 dF,o G _1 o ft'(f)] 


(15) 


and 


T = 1-4 I ^ O ^ (t) 


'0 


rp 1 oil>, 1 oFioG 1 oi/j, 1 (t) 


ft ^ i ul ft 

I- 1 / 


rp 1 °i/> 4 > 1 oF j oG 1 oip (j> 1 (t) 
y p- 1 oDj k oi/j- 1 (t) 

p' ° P" 1 - V ^ 1 °F k o Fr 1 o^o p(z)} cb 


p'°p 1 {Vv 1 W - ° o Fj 1 O Vv o p(z 

(iiH - 


-4 [ ' i 2 ^ 1 ° Fk ° G" 1 - (D fci + D w )) o 

Jo 

where 

G(z) = p {p- 1 o y- 1 o J^O) + p” 1 o y- 1 o F, (z)} + if- 1 o F fc (z) 


Proof. See Appendix O 


(16) 

□ 


Remark 3.4. The Gumbel case 

In the setting of Remark\3J\ we get 


*• ■ — * 


£(*) ^^^(-iogt ) 1 7 : 


■'(^)vl°gi) e 


A,- 


+ 


^ r>? (G-H(.-toSt) + ))* 


-logft - ^ 2 ^ 


1-4 


ofo+ 


- ( — + — ) (- log t)+ + 2A fc G- 1 ((- log tf) 

l-^jk 
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and 


T = 1 ~ 0 J o t(-logt) 1 -*- 


fA/(G- 1 ((-log^))f 

7^r) ? (-logt) s 
rjfc / 


log t) 


+ - ] dz+ 

A? 


l- 


r-Af («?-*((-log *)*))$ 


** \ Z 


'(-log t) e 


1-1 


(-log*)*-^M) dz 

A,; 


4 


7 / t(-bgt) 1 -^ 2A fc G- 1 ((-log^)^(^- + ^- 

0 \ xfJ'ik l^jk 


dt ~b 

Afc , A k 


log t)^ ) dt 


where 


G(z) — ( A?+aH V + A fc z. 


4 An application to the European banking sec¬ 
tor 


In this section we apply the model to the issue of evaluating sys¬ 
temic risk and contagion in a set of European banking systems. Un¬ 
til now, the European banking system has been segmented at the 
national level, and only after November 4th 2014 it is unified un¬ 
der a common European regulation and supervision setting (the so 
called SSM, Single Supervisory Mechanism, see, for example, Ferran 
and Babis, 2013). It is then important to recognize the relevance 
of systemic risks and contagion at the national level, and address 
the issue whether they co-move at the cross-country level. The task 
is to check whether the exchangeable contagion model may provide 
a good representation of the data. Of course, here our interest is 
mainly in the illustration of the estimation technique and how it can 
provide a guide for the specification of the model. 


4.1 Data 

We apply the model described above to a sample of 35 banks rep¬ 
resentative of 8 countries of the Euro area. The sample used is the 
same as in Baglioni and Cherubini (2013). While we refer the reader 
to that paper for an in-depth description of the data set, here we 
simply mention that the sample consists of those major European 
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banks that were subject to the stress test exercise in 2012 and for 
which a time series of CDS quotes was available on Datastream. The 
sample consists of daily data of CDS quotes, ranging from January 
2007 to end of August 2012, with the exception of Greece for which 
the sample begins on September 21st 2009, Portugal and Spain, for 
which the sample starts in January and February 2008, respectively. 
The survival probabilities were extracted from the 5 year CDS quote 
using what is called the ” simple rule”, that is assuming a flat default 
intensity, which is consistent with the model described in Section [2711 
Moreover, since it is well known that data extracted from market 
prices embed a risk premium, that is are computed under the risk 
neutral measure, we changed the default probabilities by applying 
the Sharpe ratio, according to the technique used by the Moody’s 
rating agency (see Dwyer et al., 2010). 

Future research could investigate further the marginal structure 
including more sophisticated technologies to ’’bootstrap” (in the fi¬ 
nancial literature meaning of the term) the term structure of default 
intensities (Hull and White, 2000). 

4.2 Estimation and results 

The estimation procedure applied to the data was the same de¬ 
scribed in Section [2731 Namely, we computed pairwise Kendall’s tau 
values for the survival probabilities of all the banks in the same coun¬ 
try. Then, for each country we estimated the a*, parameters and the 
6 parameter, minimizing the distance between theoretical and sam¬ 
ple Kandall tau’s. In our specific application, we used the quadratic 
distance. Due to the presence of local minima, that arose in pre¬ 
liminary work, the analysis was finally carried out using a standard 
global optimization technique, namely simulated annealing. 

Our estimation strategy consisted of three steps. 

• We first estimated the model on the whole sample for each 
country. 

• Then, for each country we used the model to estimate the inten¬ 
sity of the systemic shock and we computed the Kendall’s tau 
between the survival probability of each bank and the systemic 
shock. We verified for which countries the model specification 
is consistent with the data. 
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• For the countries where the model specification was considered 
consistent, we provided an analysis of the stability of parame¬ 
ters, by repeating the estimate in a sequence of rolling windows. 

In Table [T| we report the results for the estimates carried out over 
the whole sample. For each country, we report: i) the parameters 
for each bank; ii) the average a for the country; iii) the contagion 
parameter 9 for the country. It is worth mentioning that the average 
parameter a is computed as the harmonic mean of the a*. for each 
country. 

Before discussing the parameters estimated, we use them to pro¬ 
vide a visual check of the specification of the model. For each coun¬ 
try we provided a diagram in which on the horizontal axis we re¬ 
ported the and on the vertical one the Kendall’s tau value. In 
the diagram we plot the estimated Kendall’s tau statistics between 
the survival function of a systemic shock, estimated from the inten¬ 
sity in equation (|9]h and the survival probabilities of the banks in 
the sample. In the plot we also reported the straight line on which 
the Kendall’s tau values should lie if the model is well specified, 
according to equation (|8l). 

Figure [H shows that the model provides a good specification for 
Portugal, France, Spain and the Netherlands. In particular, the 
specification looks very good for Portugal and France. Spain and 
the Netherlands provide an interesting insight on the model. In both 
these countries there are banks whose value of ag. is very close to 
zero. They are Abn Amro in the Netherlands and Santander and 
BBVA in Spain. These three banks are not affected by the systemic 
shock, and they are very close to the intercept of the line. Never¬ 
theless, they are linked by a substantial degree of dependence to the 
systemic risk factor, meaning that an increase of their probability of 
default may be associated to a higher probability of a country-wide 
shock. 

For the other countries (see Figure [2]), that is Italy, Greece, UK 
and Germany, the model does not seem to fit the data well. In 
particular, Italy, Greece and UK could have a chance of a better 
fit if a non-exchangeable model were used. For Germany, instead, 
the model appears completely wrong, since almost all the Kendall’s 
tau’s lie below the diagonal in a region, and are not even consistent 
with the pure systemic risk specification. So, in this case, it seems 
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GERMANY 

SPAIN 

Bank 

a 

Bank 

a 

HSH 

0.288926 

PASTOR 

0.295905 

WEST LB 

0.783299 

BINTEL 

0.397199 

POSTBANK 

0.782422 

SABADEL 

0.699246 

DZ BANK 

0.637563 

POPULA 

0.842172 

BAYERN LB 

0.885879 

CAJA MADRID 

0.059619 

COMMERZ 

0.869757 

BBVA 

6.99-10 -7 

DB 

0.752793 

SANTANDER 

8.44-10~ Y 

0=1 

<5=0.625488 

0=5.803539 

<5=0.197725 

ITALY 

NETHERLANDS 

Bank 

a 

Bank 

a 

UBI 

0.260332 

SNS 

0.705556 

M-PASCHI 

0.260332 

ABN AMRO 

1.279263 

INTESA 

1.0 

RABOBANK 

0.824459 

UNICREDIT 

0.583622 

ING 

0.550825 

9=5 .589623 

<5=0.449591 

0=1.279262 

<5=0.674809 

FRANCE 

GRECE 

Bank 

a 

Bank 

a 

SOC GEN 

0.690164 

ALPHA 

0.987782 

CA 

0.469473 

EFG 

0.199412 

BNP 

0.502889 

NBG 

0.726229 

0=6.060185 

<5=0.5388 

0=3.311887 

<5=0.405181 

PORTUGAL 

UK 

Bank 

a 

Bank 

a 

LLOYDS 

0.930215 

ESP. SANTO 

0.791219 

BARCLAYS 

2.88-10 U7 

BCP 

0.695637 

HSBC 

4.08-10 iS 

CAIXA GERAL 

0.215709 

RBS 

0.769870 

0=5.369255 

<5=0.408871 

0=5.369255 

<5=0.842481 


Table 1: Parameters’ values for different banks and countries 
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PORTUGAL 


FRANCE 



Figure 1: Kendall’s tau between default times and systemic shock 


that either a model of bivariate relationships without any systemic 
shock could be preferable, or that there can be some other systemic 
risk factor missing. Actually, the discussion appeared in newspapers 
and magazines during the year of the stress testing analysis would 
suggest that the German banking system is exposed to two key risk 
factors: the first is the exposure to the so-called ”toxic assets”, com¬ 
ing from the US subprime crisis, the second, less known, is exposure 
to a specific sector of obligors, namely those linked to the shipping 
business. 

So, our estimation strategy proved able to discriminate cases in 
which the model provides a good fit to the data from cases in which 
it does not. For the four cases in which the model seems to work 
for the entire sample, we now provide an analysis of the stability of 
parameters, and in particular of the contagion parameter 9 across 
the sample. We replied the estimation using rolling windows of sev¬ 
eral lengths, even though here in order to save space we only report 
the one based on one year of daily data. An alternative more so¬ 
phisticated approach would be to consider time varying parameters 
with an estimate performed on GARCH filtered residuals. This on 
one side could be more accurate, while on the other side it would be 
inconsistent with the flat intensity assumption, calling for a proper 
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ITALY 


GREECE 



specification of a double stochastic model for each marginal inten¬ 
sity curve. 

We also performed the analysis first letting all the parameters change 
through time, and then assuming the ctk fixed across the sample, al¬ 
lowing only the contagion parameter 9 to change. The reason for 
the latter choice is twofold. First, since the sensitivity of each bank 
to shocks mostly depends on its balance sheet, it is reasonable to as¬ 
sume that the parameters a*, remain quite stable across the sample. 
Second, it was interesting to check whether the estimation of the 
contagion parameter, that is the main target of our research, was 
affected by changes in the bank specific parameters. We found that 
the dynamics of the contagion parameter is almost indistinguishable 
in the two cases. 

In Figure [3] we report the results of the analysis for the four cases in 
which the model works. The question we have in mind is whether 
the contagion parameters increased in the two crucial periods of the 
crisis. The first was in the first quarter of 2009, when the Lehman 
crisis of September 2008 propagated to Europe. The second is the 
sovereign debt crisis triggered by Greece in 2010 and then spread 
to the other countries of Southern Europe. Figure [3] confirms an 
interesting co-movement behavior of the contagion parameters for 
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Spain, Portugal, France and Greece. Intuitively, when financial cri¬ 
sis spread in the international environment, the relevance of conta¬ 
gion from the banks within each country is increasing. Differently 
from this evidence, however, in the last part of the sample, char¬ 
acterized by the Italian sovereign crisis, only contagion within the 
French banking system seems to markedly increase, while in the 
other countries it remains stable or decrease. This could be con¬ 
sistent with the greater involvement of the French banking system 
with the Italian one. In fact, in a previous version of this work, in 
which the dynamic analysis had been carried out for Italy as well, 
the contagion parameter for the Italian market was in that period 
almost indistinguishable from that of the French one. 



Time 


Figure 3: The dynamics of the contagion parameter. 


5 Conclusions and future extensions 

In this paper we presented a model that includes both systemic risk 
and contagion. Systemic risk is represented by the presence of a 
shock that brings about the default of all the elements in a clus¬ 
ter. Contagion is represented by the links between the idiosyncratic 
shocks specific to each component and the systemic shock. 
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On theoretical grounds, the analysis can be carried out assuming 
whatever dependence structure among the non observed components 
representing the shocks. On empirical grounds, here we provide an 
estimation procedure for a model in which the dependence structure 
of the unobserved components is Archimedean and exchangeable. 
Moreover, we provide a technique to verify whether the specifica¬ 
tion proposed fits the data. We also show that, including further 
restrictions may transform the copula model in a new full-fledged 
multivariate model with exponential marginal distributions. 

Given a panel data of observations of marginal intensities, the es¬ 
timation of the model is carried out on the set of bivariate depen¬ 
dence statistics. Based on estimates, one can extract the time se¬ 
ries of the systemic shock, estimate the Kendall’s tau’s of the ob¬ 
served marginals and the systemic shock, and verify whether they 
are aligned on a straight line, as predicted by the model. We apply 
this technique to a set of European banks of 8 countries, assuming 
a systemic shock at the country level, and we found that our model 
turns out to be well specified for 4 countries: Spain, Portugal, France 
and The Netherlands. For these countries, we also report an analy¬ 
sis of the dynamics of the contagion parameter, providing empirical 
evidence of co-movement in periods of international crisis. 

Of course, the next step of this line of research would call for estima¬ 
tion of the more general, non-exchangeable setting, that has been 
also formalized in this paper. More precisely, we see three main 
promising fields of development 

• Estimating the dependence structure of the unobserved compo¬ 
nents directly. Most likely, this would involve the application of 
Simulated Maximum Likelihood (SML) or similar techniques, 
in which one tries to estimate the parameters by simulating 
data as close as possible to the observed ones. Doing this can 
be very easy or very complex, depending on the degree of gener¬ 
ality that one is willing to accept. As the simplest case, assume 
one could consistently estimate the parameters in our model. 
In this case, it would suffice to estimate the systemic and the 
idiosyncratic components from the data and study the depen¬ 
dence analysis on those. Exploiting the invariance property of 
copulas, one could directly obtain a consistent estimate of the 
contagion parameters. As the most complex case, assume that 
the dependence structure of the unobserved components must 
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be handled in full generality. In that case, the concept itself 
of the a*, parameters would be lost, since there is no guaran¬ 
tee that the same proportionality between the systemic shock 
intensity and the marginal intensity is maintained through the 
sample. 

• In a similar line of research, one could also decide whether to 
focus on the full specification of the model, or only in the re¬ 
lationship between the systemic shock and the marginals. In 
the latter case, the dependence structure among idiosyncratic 
components would play the role of nuisance parameters. For 
example, our findings of exchangeable contagion could be con¬ 
sistent with a dependence structure in which some degree of 
non-exchangeability is present, but it is limited to the idiosyn¬ 
cratic shock dependence. Within this framework, estimation 
techniques such as those envisaged above could be used to de¬ 
vise formal tests of the weaker concept of exchangeability dis¬ 
cussed in the extension of our model, in which only the pairwise 
dependence between the idiosyncratic shock and the systemic 
ones are required to have the same copula. 

• Finally, on a different line of research, one could use the estima¬ 
tion procedure applied in this paper as an exploratory tool to 
identify clusters of components that may constitute the same 
’’exchangeable systemic contagion cluster”. This could be done 
evaluating the dependence between new element and the sys¬ 
temic shock representing a cluster. Or it can be obtained by 
measuring the dependence between the systemic shocks of dif¬ 
ferent clusters to evaluate if some of them can be merged in a 
single one. 

As for our specific application to the banking system, of course, the 
main challenge would be to extend the analysis to the new unified 
European banking system, represented by the 130 banks that are 
since now on under the supervision of the European Central Bank. 

6 Appendix 

Proof of Proposition 12.11 
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Proof. By (J2J), t k = 
C( u) = 0 ^ _1 (F 0 
Let 


H 0 1 of; 1 (u k ). Hence 
(max o 0~ 1 (ufc)}) ) + o F k o H~l o ^ _1 (w fc ) 


Aj = <j u G [0,1] : max{ff 0;! o 0 0q)} = tf 0J o 0 ( Uj ) 


then 


<7(u)l A ,(u) =0 [0 1 (F 0 (H 0 joip \uj))) + ^2fj 1 o F k o H 0 1 o 0 \ Uk ) 


k =1 


1 (u j ) + ^ 1 oF k° H 0 l o 0 ^Ufc) 

fe=l,fe^ i 7 


□ 


Proof of Theorem 13.11 

Proof. In the sequel we set <9i(7(w, u) = J^C(u,v) and d 2 C(u,v) = 
&C{u,v). 

We want to compute the C-measure of the set 

S t = {(w, v) G [0, l] 2 : C(u, v) < t}. 

Notice that the level curve C(u,v) = t intersects the graph of the 
function v = h(u ) in a unique point that we denote with (u t ,v t ). 
Hence S t can be decomposed as S t = Rt + Ri,t + R 2 ,t where R t = 
[0,tq] x [0,v t }, R lit = {(«,u) : v G (v t ,l],C(u,v) < t} and R u = 
{(it, v) : u G (u t , 1], C(u, v ) < t}. 

Clearly, the (7-measure of R t is t. In order to compute the C- 
measure of R.\ t and R 2t , we compute u t and v t . Since (u t , v t ) satishes 
% + D ki ° iff X {v t )) = t and v t = h(u t ), we get 

iff 1 °^e° Df 1 o D ik o iff l fvf) + D ki o = iff 1 ft) 

from which 

v t = 0^oG^ 1 o0- 1 (t) 
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and 


u t = i> e o D i3 o D ik o G o i/) (j> (t). 

Let us start with Ri t . Notice that here, C(u, v) < t is equivalent to 
u < Fi(t, v ) where F^t, v ) = ^ (^(t) ~ D ki o^/(u)). Hence 

P(i2 lit )= f F(U <F 1 (t,v)\V = v)dv = 

J Vt 

= f d 2 C(F 1 (t, v),v)dv = 

J Vt 

r 1 d 

= J ’ -j-DkiOip^^dv = 

= ° ■ D ki o 

Let us now consider R 2t . Notice that here, the inequality C(u,v ) < 
t, is equivalent to u < F 2 (t,v) where 

F 2 (t, v)=^ e o Djl ( p- 1 (^\t) - D ki o ^\v)) - D ik o ^ l (v)) . 

But 

R 2 .t = {(u, v) : u t < u < 1, t < v, C(u, v ) < £}U{(w, v) : u t < u < 1, v < t} 
and 

P (ut < U < 1, V <t) — t — C(ut, t ) = 

= P ([/ < u t , t < V < v t ). 

Hence 

/ Vt 

F(U < F 2 (t,v)\V = v)dv = 

/ Vt 

d 2 C(F 2 (t,v),v)dv = 

pVt 7 7 

= / ^ ° ° P _1 (%\t) ~ D ki ° v )) D ik o v ) + — D ki c 


dv 


"Dik 0 ^^ 1 (vt) 


P °P {% (t) ~ D ki ° D ik ( z))dz+ 

D ik °4> 0 1 (t) 


+D ki o ip^ivt) - D ki o t )} . 

From P (S t ) = t + P(i?i it ) + P(i? 2 ,t) we get (TT5|1 


tl>4,\v)}dv 
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As a consequence, the Kendall’s tan is 


r = 3 — 4 f K(t)dt = 

Jo 


= 3-4/ {t-^o^^t)- 
Jo 

r-DikoG^oip- 1 ^) 

D ki oii~ l (t)- / 

J DikOip^tt) 


p'°P V/0) ~ D ki° D ik{z))dz 


dt 


= 1 + 4 / %°% (*)■ D ki ° % ( t)dt- 
Jo 

rD ik oG- 1 oi,- 1 (t) 


4 / 


' Dikoip^it) 


p'° P 


□ 


Proof of Theorem 13.21 

Proof. In the sequel we set d\C(u,v) = y -C(u,v) and d 2 C(u,v) = 
mC(u,v). 

The proof is similar to the one of Theorem 13.11 

Again we decompose the set S t = {( u,v ) G [0, l] 2 : C(u,v ) < t} 
as S t = Rt+Ri,t+R 2 ,t where, if (u t , v t ) is the intersection point of the 
curves C{u, v) = t and v = h(u), R t = [0, u t \ x [0, v t \, R\ ]t = {(w, v) : 
u G (ut, 1], C(u, v ) < t} and R 2>t = {(u,v) : v G (v t ,l],C(u,v) < t}. 
Clearly, the C'-measure of R t is t. In order to compute the C-measure 
of Ri t and R 2} t, we compute u t and v t . Since 

Vty [p (d” 1 ° Afc ° Vt) + P~ l ° D jk ° V^V)) + D kj o = t 

and v t = h(u t ), we get 

ft = V ° H 0jj o (T _1 o ^ 1 (t) 

and 

= Vty o #o,i ° G -1 o 
where G* is given by (jT6T) . 

Let us start with Notice that here, C*(u, v) < t, is equivalent to 
u < Fi(t,v) where 

Fi(t, v) = ^oDr^op (p- 1 (^(t) - D kj o ip- 1 (n)) - p" 1 o D jk o . 
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By similar arguments as those used in the proof of Theorem 13.11 we 
have 


rvt 


P (R ltt )= V(U<F 1 (t,v)\V = v)dv = 




d 2 C(F 1 (t,v),v)dv = 




d 


P'°P 1 (%\t) - DkjO^iv)) ^P 1 °D jk oiP (j) 1 (v)+ 


d 


+ di’ Dkj ° ^ ^ f dv = 


P 1 °Dj k oil>' j> 1 (vt) 


+ D kj ° ~ D kj O . 

Substituting v t we get 

P(^i,t) = 

= 


P°P oFkoFj o^o p( z ))dz+ 

P~ 1 oD jk op (j) 1 (t) 


p 1 oil> <t> 1 oF j oG 


■ lu ^ J / p’ O p ^/(f) - o F k O F j 1 O Vv o pO))dz+ 

p - 1 oD jk op r/j 1 (t) 

°Fk o G _1 o - D kj o ^ 1 (t)} ■ 

With similar computations we get 
P(i22,t) = 

| />p _1 


CP 1 opPoF i oG 1 tp 6 1 {t) 


'• J / _ P'°P 1 (^ 1 W -^4 > 1 °FkoF i 1 o ^ o p(z))dz+ 

1 oD ik oip <j> 1 (t) 

°F k ° G' 1 o ^(i) - ° • 

From P(5' t ) = t + P(i?i it ) + P(i? 2 ,t) we get (1T5]) . 
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As a consequence, the Kendall’s tau is 


r = 3 — 4 / JC(t)dt = 

Jo 

= 1-4: [ 

Jo 

rp~ 1 °i>^ 1 °F'joG- 1 oi,- 1 (t) 


P 1 o*l>. 1 oF i oG 


J p~ 1 oD jk o^ <j> 1 (t) 

0 1 oFioG 

’P~ 1 oD ik o^~ 1 (t) 

r l 


P °P 1 {V’/W ~ ° F k ° Fj 1 ° Vv ° p{z)} dz+ 


f> ° P 1 {% 1 ( t ) ~ 1 ° F k ° 1 ° VV ° P(-)} ^ 


(itH - 


-4/ ip# o ^ (t) ■ ( 2^ oF k oG - (D kj + D ki )) (t)dt 


□ 
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